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For one-dimensional vibrating cavity systems appearing in the standard illustration of the dy- 
namical Casimir effect, we propose an approach to the construction of exact closed-form solutions. 
As new results, we obtain solutions that are given for arbitrary frequencies, amplitudes and time 
regions. In a broad range of parameters, a vibrating cavity model exhibits the general property of 
exponential instability. Marginal behavior of the system manifests in a power-like growth of radiated 
energy. 
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I. INTRODUCTION 

The dynamical variety of the Casimir effect has been studied in numerous papers in the last two decades [l], 0] • 
The most favorable model for theoretical considerations refers to the original Casimir's setup [3J of a cavity composed 
of two perfectly conducting parallel plates. The dynamical modification means that the distance between plates 
changes with time. Then, we come across new fascinating phenomena in the context of quantum field theory. The 
most impressive manifestation of unusual non-classical properties of the theory of quantized fields is the effect of 
particle production "from nothing". It is referred as dynamical Casimir effect (DCE) or motion induced radiation 
(MIR). Most of theoretical papers explore so-called vibrating cavities, where oscillations of cavity walls are periodic 
in time. Attention of their authors is attracted by the instability of quantum fluctuations due to the parametric 
resonance. They hope the resonance enhancement of particle production could lead to macroscopic effects that would 
be experimentally detected. Nowadays, the resonant instability of the vacuum that is followed by an explosive particle 
production in a vibrating cavity is believed to be accessible for measurements [1, H, H 0, 111 . However, we are still 
awaiting for confirmation of successful experimental projects. 

In this paper, we consider the behavior of the electromagnetic field in a one-dimensional vibrating cavity. This 
problem was studied initially by Moore [9j . A cavity with one stationary wall and one moving wall with some prescribed 
trajectory had been elaborated there. Moore did not solve any particularly interesting cavity model, but he made a 
thorough study of a general theory of solving such models. His basic approach was generally used and developed in next 
years [To|-[7^|. From theoretical point of view, this simplified model deals with several hard and important problems. 
We mean the task of solving a wave equation with time-dependent boundary conditions, the difficulties with analytical 
description of physical systems under the parametric resonance conditions, the problem of quantization of fields in 
limited regions with moving boundaries, the squeezing of quantum states or the problems of quantum entanglement 
and decoherence. One-dimensional vibrating cavities provide the simplest theoretical laboratories to study these 
issues in quantum field theory. The methods that are worked out there can be essentially adapted for more complex 
models. Some results and ideas are endorsed as well. There are advanced approaches to include non-perfect [H, [56j 
or p artly transmit ting [72 1 cavity walls, finite temperature effects [U Hfl [47j or proceed to three-dimensional case 
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In the case of one-dimensional cavities and " scalar electrodynamics" , main achievements of numerous investigations 
were obtained either in the framework of the effective Hamiltonian approach [l|, l23l |35| or using various numerical 
approaches [24|, |48|, |5J, [66|, [67|, |69|, |70j . Analytical solutions obtained through the perturbation methods with effective 
Hamiltonians hold only for small amplitudes [H, US H3 or for particular time regimes, either short time [3(| or long 
time [20] limits. In many investigations, a frequency of cavity vibrations is assumed to match a resonance frequency. 
It is well known from classical mechanics [73j that parametric resonance occurs also for frequencies that are not finely 
tuned provided that respective amplitudes of oscillations are sufficiently large. Off resonant behavior of vibrating 
cavities is usually studied in the limit of small detuning from resonance frequencies [38l l58j . 

Our aim is to gather exact analytical and global solutions for vibrating cavities. In fact, there are few known 
solutions that can be described by closed form expressions. The first closed form exact solution to describe a cavity 
vibrating at its resonance frequency was presented by Law [2lJ . Law's solution corresponds to a cavity that oscillates 
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basically sinusoidally for small amplitudes. The frequency of oscillations is twice the lowest eigenfrequency of the 
cavity (so-called "principal resonance" [J]). Law found travelling wave packets in the energy density of the field. He 
noted " sub-Casimir" quantum fluctuations far away from the wave packets. Next, Wu et al. presented a family of 
exact analytical solutions for all resonance frequencies. In the particular case of the second resonance frequency, their 
solution is matching Law's solution. They described emerging of wave packets in the energy density, indicated sub- 
Casimir fluctuations and emphasized the absence of wave packets in the first resonant channel with the fundamental 
cavity eigenfrequency ("semi-resonance" One can be puzzled that for any one of known exact solutions, there 
appears a power-like resonance instability. Faithfully, the total energy of the field increases quadratically with time 
there. On the other hand, it is well recognized from other cavity models that the total radiated energy typically 
grows exponentially with time 0, H3, H3, IH, 0, EH • For instance, one can refer to the asymptotic formulas found by 
Dodonov et al. [l[ for cavities that undergo harmonic oscillations. It was generally argued [3[27| that an exponential 
resonant instability is typical for vibrating cavities, while a power-like behavior constitutes a critical boundary between 
stability and instability regions defined by domains of parameters 63]. In this paper, we find solutions that reveal 
exponential instability generally and exhibit a power-like law as a marginal effect. Moreover, all previously known 
exact solutions applied only to resonance frequencies. Here, we provide exact solutions that are adjustable for all 
frequencies. This paper presents a rich class of exact and closed-form solutions, in addition all formerly presented 
solutions [l|, I2H [41 1 are captured here as particular cases and examined in a more comprehensive way. 

Our paper is organized as follows. In Section II, we present our way of representing solutions to describe the 
quantum dynamics in a vibrating cavity. It relies on SL(2,R) symmetry of the algebraic structure that exists for 
the quantized scalar field in a static cavity [5(1 HH [zl] • Actually, we abandon Moore's function that is awkward in 
use. We put forward another object, called a fundamental map by us, that has remarkable analytical properties. 
These properties are collected in Section III, together with appropriate mathematical formulae for primary physical 
functions, namely the vacuum expectation of the energy density of the field inside a cavity and the total radiated 
energy. Finally, in Section IV we present a collection of exact closed-form solutions. The results are summarized in 
Section V. 



II. REPRESENTATION OF SOLUTIONS TO DESCRIBE THE QUANTUM DYNAMICS IN A 

VIBRATING CAVITY 

In the standard physical setup, we have an electromagnetic resonator of length L composed of two perfectly reflecting 
walls. Initially, the cavity is static. Then, it undergoes vibrations with a constant frequency uj. In literature, it is 
frequently assumed that the cavity length L is related with the period of oscillations T — 2tt/lo. In this paper, we 
will keep the parameters L and T independent. The parameters provide the characteristic physical length scales. The 
static cavity length L defines the magnitude of Casimir interactions. In particular, it specifies the scale of quantum 
fluctuations leading to the production of particles. The period T is the scale of parametric excitations of the system 
caused by some external force. It is very useful in numerical computations to put T — tt. The parametric resonance 
is expected when L and T are of the same order. Eventually, it depends also on an amplitude of vibrations. We are 
willing to yield a phase diagram (Arnold's diagram [73[) that exhibits stability and instability regions. 

The derivation of the simplest mathematical model leads to the quantization of free scalar field A(x,t) with Dirichlct 
boundary conditions imposed at the boundary walls x = and x — L(t). The trajectory of the oscillating wall is 
periodic: L(t + T) = L(t). It is important to assume that L(t) > (the cavity never collapses) and \L(t)\ < v max < 1 
(the wall velocity does not come near the speed of light). Moreover, we impose L(t) = L for t < (the cavity is static 
in the past, this condition is important for the quantization). The construction of the basic set of solutions for this 
problem was given in Q: 

i 

A]y(t,x) — = [exp (— iu>pjR(t + x)) — exp (— ioj]^R{t — x))] . (1) 
\JA_itN 

The cavity eigenfrequencies lo^ = Nir/L are called resonance frequencies. We expect that the parametric resonance 
occurs at these frequencies for any amplitudes. However, the instability of the system may appear also for other 
frequencies provided that the amplitude of oscillations is sufficiently large. Usually, it is a hard task to get the picture 
of the asymptotic behavior of the system for any frequencies and amplitudes. Our knowledge of the system comes to 
us through the Moore's function R given by the following equation: 

R(t + L(t)) - R(t - L(t))=2L . (2) 

Usually, Moore's function R is defined as a dimensionless function (phase function). In this paper, we will prefer to 
define this function in dimensions of length. There is no general theory of solving Eq. @ . Before we present a big set 
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of exact solutions of the above problem, it is worth to recall some useful symmetry of the static cavity system [50|, l2J] . 
In the static region for t < 0, the quantized theory is invariant under the conformal transformations: 

t ± X — > Rmin{t ± x) , (3) 

with the functions R m in defined by: 

Rmin{T) = — arctan (cr (tan ^i^)) , (4) 
uj\ V 2 / 

where cr(r) = (ylr + B)/(Ct + D) is any homography and io\ is the lowest resonance frequency. Subsequent branches 
of multivalued function arctan should be always chosen and linked together in such a way that a resulted function 
Rmin is continuous. It is described here a well-known SL(2, R) symmetry of free scalar fields quantized on a strip 
[74j . Surprisingly, this symmetry is rarely exploited in numerous papers on physical models of the quantum field in 
a one-dimensional cavity. In particular, the symmetry helps to solve the puzzling problem why there is no resonant 
behavior of the system for the fundamental resonance frequency io\ = tt/L. 

In this paper, we will be searching for exact solutions of Eq.Q in the following form: 

i?(r) = — arctan ^A„( T -) (tan + shift. (5) 

In order to obtain closed-form solutions, we assume the range of arctan to be [— tt/2, tt/2] (principal branch) and 
appropriate shifts will be explicitly specified throughout. For instance, the linear Moore's function, which describes a 
static cavity, should be represented as: 

. . 4t 2 lot lot 3, 2ir 

RstaUc(T) = t = — arctan (tan — ) + [- -J — , (6) 

to to I Zir I u> 

where we have used the standard notation for the floor function. The construction of the representation Eq.([5]) is 
tied up with the well-known idea from classical mechanics (73|. To explore the dynamics of periodic systems with 
parametric resonance, it is a handy way to deal with mappings for single periods. Here, we need a set of maps A n 
numerated by the number n. Fortunately, the maps are not independent. We prove that it is enough to specify only 
the first map Ai. Henceforth, a function Ai(y) is going to be called a fundamental map throughout this paper. This 
map defines the auxiliary function /: 

/(r) = — arctan fAi (tan — )] + shift, (7) 
lu V 2 / 

which is a solution of a simpler problem than Eq.Q (see equations for billiard functions in J63|): 

f(t + L(t))=t-L(t) . (8) 

Since the cavity is static in the past, we have always that /(r) = r — 2L for r < L. The subject is also simplified due 
to the fact that the auxiliary function / fulfils the periodicity condition: 

f(r + T) = f(r)+T . (9) 

In general, the Moore's function R(t) is not subject to any periodic conditions. The reason lies in the lack of periodicity 
of the index n(r), that assigns a map to a particular point. 

It is straightforward to prove that a fundamental map Ai(w) designates unambiguously a Moore's function R(t). 
The solution of Moore's equation (J2]) can be build according to the formula: 

R( T ) = /°«W( T ) + 2i[n(r) - 1] . (10) 

Looking at the representation Eq.©, one can check easily: A„ = (Ai) 0,i . Throughout this paper, we use (Ai) on to 
note n-fold composition Ai o Ai o ... o Ai. It remains only to describe the step function n(r) that appears in Eq.© 
and Eq. (TlT)|) . As the function /(r) is increasing, the region for r > L can be covered by intervals [i n _i,i n ), where 
L n = (/ -1 ) 0,i (L). The map number n(r) equals n if the point r lies inside [X n -i, L n ). Map markers L n will be called 
milestones throughout this papers. If r G [L n -i, L n ), then /(t) 6 [L n —2, L n -x). Thus, it is easy to find the following 
recurrence relation, which is also very convenient for numerical purposes: 

" (T) = {l + n(/(r)) T tVl ( H ) 
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In order to provide a glimpse to details of future calculations with the representation Eq. ([5]) or Eq. ([7]) , we take a look 
at the solution given by Eq.(|4"5|). This solution will be discussed later, but we glance over the borders of intervals 
for corresponding functions / and R there. They are depicted in FigQ] Performing appropriate calculations, one 
should take into account that all variables and mappings are valid only in defined domains. Typically, the region 
suitable for calculations of Bogoliubov coefficients or total radiated energies is covered by two subsequent maps A„. 
It makes evaluations of integrations and derivations of formulas more complex. This is the price we have paid for 
the replacement of Moore's equation by simpler relation Eq.®. However, we will convince ourselves that this way 
is effective as a method for obtaining analytical results. Later, the details of calculations will be always skipped, so 
that the pattern in Fig.l is the only commentary on practical calculations. 
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Figure 1: The borders of different change intervals for the Moore's function R defined by Eq.Q with Eg. ([48 



It is difficult to derive the function R(t) from Eq.([2]) for some prescribed trajectory L(t). A great number of 
numerical approaches and approximate solutions were presented in other papers, but only few exact solutions are 
known. One way to obtain exact solutions is to specify the function /(t), and then the trajectory L(t) can be given 
in a parametric form: 



* = [t + /(t)]/2 
L(t) = [r-/(r)]/2 



(12) 



The prescribed function /(r) represents an admissible physical trajectory provide that it fulfils the following require- 
ments 1481: 



(0 
(ii) 
(Hi) /(t) 



l—v„ 



t - 2L for t < L 
< fa) < 



(13) 



< T 



In this paper, we will be exploiting the representation Elq. (|5j) to describe solutions of equations for the electromagnetic 
field in an oscillating one-dimensional cavity. Before we start with the construction of solutions, we describe general 
properties of fundamental maps Ai extracted from proper solutions. 
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III. GENERAL PROPERTIES OF FUNDAMENTAL MAPS Aj 



Knowledge of Moore's function enables us to draw out all information about the vibrating cavity system. The most 
important object to calculate is the vacuum expectation value of the energy density: 

(T 00 (t,x)) = Q{t + x) + e {t-x) . (14) 

Using appropriate formulas given in [64| and our representation Eq.([5]), we can easily calculate: 

,2 ,.,2_,.,: 



q(v - 



48tt 48tt 



1 + 1,2 A' f \Y 



1 + 



—(l + v yS[A n{T) }(v) , (15) 



where v = tan (lot/2) and S[A n ^](v) denotes the Schwartz derivative of A„( T ) with respect to v. The total quantum 
energy radiated from the cavity can be calculated from: 



f L(t) r t+L(t) 2 f dv 

E{t)= / dx(T 00 (t,x)) = / drg(T) = - / — - g(v) . (16) 

Jo Jt-L(t) U J l + V z 



The most useful is the last formula which enables us to calculate the total energy by integration with respect to v. 
However, we should remember from the comment in the previous section on the pattern in Fig[T]that the replacement 
v = tan (lot/2) is valid only for a single period of cavity motion. The interval of integration [t — L(t),t + L(t)} is to 
be divided into parts representing separate periods of motion. The map number ti(t) may change at most once per 
period. 

Let us remind the relation A n = (Ai)°™, and we need only to specify the fundamental map Ai. The knowledge of 
this map makes it possible to predict the evolution of the system and describe the resonance behavior. Henceforth, 
our exploration of a quantum field theory system is quite similar to examination of classical mechanics models under 
the parametric resonance [73| . We need only to analyze the asymptotic behavior of iterations of the mapping ruled 
by Ai, which is known from the first period of motion. 

We are going to make a list of general properties of fundamental maps Ai. First, we include information that the 
cavity is assumed to be static in the past, i.e. for times t < 0. It follows that a solution for t > is to be sewed 
together with a static one at t = 0. In our context, there is no need to demand that the sewing is perfectly smooth. 
For instance, we can accept that a force which causes cavity motion may be suddenly switched on. Such a solution 
may lead to some Dirac delta terms in its function for energy density, but from physical point of view the solution 
is acceptable and useful for applications, so that it is definitely worth saving them. Henceforth, let us propose some 
minimal set of requirements for sewing. We put forward three sewing conditions at the initial time t — 0. The 
trajectory of the cavity wall and its velocity should be continuous: L(t = 0) = L and L(t = 0) = 0. Moreover, there 
should be no sudden local growth of energy: (Tbo(i = 0, x)) = — tt/(24L 2 ), i.e. the local energy density matches the 
Casimir energy density of vacuum fluctuations at the initial time. It is now straightforward to gather a full set of 
initial conditions for the fundamental map Ai(v): 

= ~ V ° loL f L\ fn lo\ 

A 1 (v ) = 1 ; v = tan — = tan tt— = tan - — . (17) 

S[Ai](«o) = 2 V TJ V2^i/ 

The last condition implies that the construction of a fundamental map is yet a non- linear problem. Next, we impose 
the requirement that the velocity of cavity wall should never exceed v max . The maximal velocity is a parameter of 
the cavity model and the only limitation is that v max < 1< From Eq. (|13p (ii) we obtain: 

) - Vmax < , l \t M v) < ] + Vmax . (18) 
1 + v max ~ 1 + A\(v) iyJ ~l-v max 

This is a strong constraint on possible maps. One immediate consequence is that our function is increasing: A[(v) > 0. 
We can also learn about its singularities from Eq. (fl8|) . If the function Ai(v) is singular at some v s : 

lim Ai(u) = ±oo , (19) 

v — >v s =F0 

then it is easy to prove that the following limit is finite and different from zero: 

lim ^4 = - lim 7 \ a 7 ; ■ (2°) 

v—m, A±(v) v—>v B (v — v s )Ai(v) 
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Henceforth, it follows the function Ai(u) may have only poles of order one: 

M») = 7 w ^ 7 \ > (21) 

(w - Vi)(v - v 2 )—{v - v s ) 

where the numerator h(v) is an analytical function. Taking Eq. (fl"8|) together with Eq. (|2"Tj) , we note that for large 
values of v the function h(v) shows the following asymptotic: 

h(v) ~ \v\ k , k G {s-1, s, s + 1} . (22) 

Finally we look at the representation Eq.® and the periodicity condition Eq.©. We conclude that the number of 
singularities s in the map Ai for the representation Eq.® is at most one. Actually, we could replace lu with suj in the 
representation Eq. 10 and allow for more complex form defined by Eq. (|21|) . The same performance as that in Section 
IV might give new exact closed form solutions, but we will not examine this idea here. 

In general, for large arguments either the function Ai is unbounded or it takes a finite limit. Therefore, the 
respective continuity condition corresponds to one of two choices: 

Ai(±oo) = ±oo or Ai(— oo) = Ai(+oo) = finite value . (23) 

Let us summarize this section. The basic set of solutions Eq. ([T]) for a quantum cavity system can be fully specified 
by Moore's function Eq.®. In turn, this function is to be reconstructed from the fundamental map Ax- The 
fundamental map is associated with the first period of motion. Some basic physical requirements lead to strong 
mathematical conditions on the application of function Ai to cavity models which are admissible from physical point 
of view. This includes suitable sewing conditions Eq. (fT7|) at some distinguished point vo, the inequalities Eq. (fT8|) 
introduced by a limitation v max on a cavity wall velocity and continuity condition Eq. (l23D . Moreover, the function 
Ai may have at most one singularity (only a simple pole) and it behaves for large arguments according to Eq. ([22|) (s 
is a number of singularities, i.e. or 1 here). 

For some given fundamental map Ai that fulfils all required mathematical conditions, it may be still difficult to 
derive trajectory L{t) from Eq.® and Eq. (|12[) or map ranges L n and index function n(r) from Eq. ljlip . However, 
it is possible for rational functions. In the following section, we will discuss such solutions. They form a big and 
interesting family of exactly solvable cavity models. In particular, they include all examples of exact closed form 
solutions on vibrating cavities known from other papers [2l| [44| . 



IV. EXACT CLOSED FORM ANALYTICAL SOLUTIONS 



We use the considerations of the previous sections to find a family of exactly solvable quantum models of vibrating 
cavities. The static cavity length L is fixed and it characterizes a physical scale. According to the naive understanding 
of parametric resonance, the frequency of vibrations u> should be close to one of the resonance frequencies ojn- It 
means that L is close to NT/2. However, it should be confirmed in a specific cavity model whether this naive criterion 
of resonance is justified. Moreover, it turns out there is a more subtle situation when L is an odd multiplicity of T/2, 
i.e. the parameter vq in Eq. (|17p is infinite. Such cases should be analyzed in our treatment separately. 



A. Linear fundamental maps Ai 

We begin by addressing the case when a fundamental map Ai is a polynomial. The condition Eq. (TrS]) that velocities 
are not approaching the speed of light is very restrictive here. It allows only for a linear function. First, we will examine 
the case when vq is finite. 



1. Finite values of vq 



Our method of proceeding follows closely on the formalism presented in the previous sections. Inserting a linear 
function into conditions (|17p . we pick out: 

Ai(v) = v-2v ; v ee tan . (24) 



7 



It is easy to verify that the above function fulfils all physical requirements Eq. (fTTf , Eq. (fTSf and Eq. (|23f . Let us define 
the natural number M and the angle parameter 9 by: 



L=(M + -)T ; M = l,2,3,... , \9\ < - 

TT Z 



(25) 



The parameter M can be interpreted as the order of the resonance. We will go through this subsection and see that 
the parameter M is better to characterize the resonance channel than N. The auxiliary function / for r > L from 
Eq.([7]) and its inverse function f^ 1 for r > — L yield: 



/(r) = 5arctan(tan^-2i;o) + ([£ + ±J - 2M) T , 
/ _1 (t) - £arctan(tan^ + 2i; ) + (|_? + |J + 2M) T , 
v = tan ^ = tan 9 . 



(26) 



The corresponding trajectory of the cavity wall for t > is to be evaluated from Eq. (fT2|) . Using some trigonometric 
identities, we reveal the following path: 

1 9 
L{t) = L-\ arcsin (sin 9 cos (uot)) . (27) 

For small parameters 9, the oscillations of the cavity wall are close to a sinusoidal wave (see Figj2]). With increasing 9, 
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Figure 2: The trajectories of the cavity walls for the cavity motion Eq. (|27[) with M = 2 and 9 = 7r/4. 

they are nearer to a triangle wave. The wall oscillations take place between MT and L. The amplitude of vibrations 
is AL = 2\9\/oj, and the maximal velocity yields: 

Vmax = | Sill 1 . (28) 



(29) 



The Moore's function for r > L can be calculated from Eq.JTU] 

2 



R(t) — — arctan ^tan — 2n(r)tan0^ + shift . 
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Figure 3: Moore's function for the cavity motion Eq, (|27|l with M = 2 and 9 = 7r/4. 



The representation Eq. ([29|) can be effectively used if we are able to assign appropriate maps. It is nice that we are 
in a position to calculate the milestones L n and the map number n(r) from Eg. (|1 1|) exactly: 

L n = 1 arctan (2n + l)v + (2n + 1)MT , 

n(r) = n (r) - 1 + 9(r - L no ( T )_i) + 9(r - L no(T) ) , (30) 
no(r) = \r/(2MT) + 1/2J , 

where the Heaviside step function is defined with 0(0) = 1. The Moore's function R(t) from Eq. (f29|) for some specific 
motion of type Eq. (f2"T)l is shown in FigJ31 This function is always a small deviation from the linear function Eq. © that 
describes the static case. Disturbances caused by the cavity motion are magnified in Fig[3]for small and large function 
arguments. With increasing arguments they approach a well-known staircase shape ("devil's staircase"). However, 
the steps are hardly regular. If we look at the asymptotic behavior of the Moore's function, then we are convinced that 
it is not a good object for practical calculations, both analytical (perturbation methods) and numerical. Thus, the 
transformations for phase functions like Eq.(jn]) are necessary to get a feasible way to perform mathematical analysis 
of vibrating cavities in quantum field theory. 

The shape function for the energy density Eq. (fT5)) for the solution Eq. lpM)) reads: 



P(r) 



48tt 48tt 



1 



1 + (v- 2n(r) tanf 



(31) 



A snapshot of the energy density is displayed in FigJU In general, there are M wave packets travelling left and M 
wave packets travelling right. Their localization and their evolution can be easily derived and it is in full agreement 
with results of procedures described in [63[ and generalized in [68[ . One can successfully derive periodic optical paths 
and calculate cumulative Doppler factors, cumulative conformal anomaly contributions and other quantities. Here, we 
skip such details. Far from the narrow packets, in the so-called sub-Casimir region [2l| the energy density is constant 
and its asymptotic value is: 



2 f 29^ 2 



T^ t {T)^~— = ~[2M+-) po , for larger, (32) 

Z47T \ 7T 
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Figure 4: The energy density of the cavity Eq. (|27|l with M = 2 and 8 = n/4 defined in terms of the Casimir energy density of 
the static cavity po = 7r/(24L 2 ). 



where po = 7r/(24L 2 ) is the magnitude of Casimir energy density for a static cavity of length L. Most of the energy 
is concentrated in narrow wave packets. The heights of peaks are proportional to i 4 , and their widths shrinks like 
i~ 2 . It suggests that the total energy grows with time like t 2 . It is true, and one can calculate from Eq. fTB)) an exact 
formula. Here, we give only an asymptotic formula for large times: 

(tang)2 * 2 ' t>>l - (33) 

As usual, there is no resonant behavior for the the lowest resonance frequency. However, the resonance emerges for 
all frequencies above this threshold: u> > uj\ (or equivalently: for L > T/2). Paradoxically, the resonance appears 
here for all frequencies but resonance ones. For resonance frequencies u>n, either the cavity is static or the motion 
is singular (a triangular wave trajectory). Therefore, we should learn that the resonance frequencies are auxiliary 
objects, and real behavior of any physical cavity system depends on its individual features. A specific feature of the 
solution Eq. (|27[) is that for some fixed initial cavity length L, the resonance appear for almost all frequencies above 
some threshold. However, there is no possibility to adjust the amplitude of vibrations. There exists an exact relation 
between the amplitude and the frequency: 

AL 

The most important problem for any linear dynamical system that exhibits a parametric resonance phenomenon is 
to find stable and unstable regimes for periodically excited parameters. Usually, the parametric resonance domains 
are depending on three crucial parameters: frequency and amplitude of periodic excitation and damping coefficient. 
The relevant Fig|5] exhibits the phase diagram for the cavity model described in this section. Since the frequency 
and the amplitude are related by Eq. ([M|) . the instable solutions Eq. ([2"T]) are represented by points on the curve to 
this plot. The instability of solutions is quadratic according to Eq. (|3"3")l . It was suggested in [63| that a cavity model 
with a power-like instability appears as some boundary limit. If we extended the model Eq. ([2T|) to possess more free 
parameters, then other points in Fig[5] would represent cavity motions when the amplitude and the frequency do not 
match Eq. (j3~4")l . Below the border curve, the cavity system would be stable, and for states represented by points that 
are placed above the curve we would observe the resonance with the exponential growth of the total radiated energy. 
The last example of solution considered in this paper will justify such predictions. However, we are not able to prove 



2u>i , LU 
io 2u>i 



(34) 
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Figure 5: The phase diagram for the cavity model Eq. (|27|) : the relative amplitude of cavity oscillations versus the frequency 
as the multiplicity of the fundamental frequency. 

that the statement is generally true. Note famous Arnold's tongue structure [73| in FigO However, the tongues are 
rather broad. In the classical theories [73j], Arnold's tongue has usually a narrow knife shape. 

2. Infinite values of vq 

We get infinite values of sewing points in Eq. (|17p if the cavity oscillates at odd resonance frequencies u> — u>2M-i, 

L = (M-i)T ; M = l,2,3,... . (35) 

The conditions Eq.|T7|), Eq. (jl8p and Eq. (|23p are satisfied by a linear map with an arbitrary intercept parameterized 
by 6 (warning: 9 has different meaning that the same parameter in the previous subsection. For convenience, we have 
redefined this parameter here in such a way that numerous formulae match those of the previous section): 

Ai(v) =«-2tan0 ; |0| < | . (36) 

The trajectory of the cavity wall L(t) is the same as in Ea.(|27p. But now, the physical situation is different. In 
the previous section, we were almost free to adjust the frequency of the oscillations. If the frequency were fixed, 
the amplitude would be given by Eq. (|34[) . Here, the frequency is not arbitrary, but the amplitude of the oscillations 
2|0|/a> is adjustable. These solutions are already known and they were presented first in [4l| (they correspond to the 
solutions numbered by m = 2N — 1 using the notation of that paper) . Our auxiliary functions are slightly modified: 

/(r) = 1 arctan(tan^ + 2tan0) + (|_fr + |J - 2M + l) T , 

(37) 

/^(t) = ^arctan(tan^ - 2tan0) + + |J + 2M - l) T . 
Again v max — \ sin0|. The milestones L n and the map number n(r) are given by much simpler formulae: 

L„ - (2n+ 1)L, 

(38) 

"(r) = L^ + 3J- 
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The Moore's function is given by the following formula: 

, . 2 [LOT 

R(t) = — arctan tan 

w lu \ 2 

The profile function for the energy density is given by: 

(2M-l) 2 vr M(M — 1)tt 



■2n(-r)UmH) + [1 + IjT - 2L . 



P(r) = - 



48L 2 



12L 2 



l + (w + 2n(r) tan (9) 2 



(39) 



(40) 



Now, it is much more easier to work out the integral Eq. (|16p . Actually, it is interesting to look into an exact and 
closed form formula for the total energy produced inside the cavity: 



jp f .\ _ M(M-l)fftan 2 8 ,% 

m? 1 + 



M(Af-l)tan 2 
3(2M-1)L 2 



sign ( 



arctan 



1 l + ^/l+sin 2 uit tan 2 8-(t/L + l-2a(t)) sin wt tan 8 

tan (wt/2) i_^i +s i n 2 ut tan 2 8 + (t / L+l-2a(t)) sin art tan 8 



(t + (1 - 2a(i))L) + 



(41) 



(2M-1) 2 tt T f,\ , M(M-1)tt . M(M-l)n tan 2 8 n\\ , 

24L 2 + 6L + — 3L ~ a ( t )) + 



2AL? 

M(M-l)tan8 
3(2M-1)I 



where 



h2(t/(2L)-q(t)) 2 tan 2 9+(t/L+l-2a(t)) tan (i^(t+L(t))/2) tan 8 
l+(tan (oj(t+L(t))/2) + (t / L+2-2a(t)) tan 6>) 2 



«(*) = — - I— J 

V / or Lr)j-J 
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-2L- 



(42) 



Similarly to the solution presented in the previous subsection, the total energy of the system grows quadratically 
with time. We have extracted the leading term. However, the next to leading terms that are linear in time play an 
important role as well. They cause that the energy is not irradiated continuously but rather in sudden jumps. To 
verify that, we should take two leading terms from Eq. (|41[) and make the approximation for large values of t. As a 
result we obtain: 



E(t) = 



M(M — l)7rtan 2 9 
3(2M - \fh 



LyJ+e(0)J 



(43) 



The presence of Heaviside function means that the problem is not analytical in the parameter 9, i.e. with respect 
to the change of direction of oscillations. We see that impulses of energy growth occur every period. For small 
amplitudes, the energy is proportional to the square of the amplitude. This is in a agreement with a non-relativistic 
limit of small velocities. It is amusing to consider a quasi-classical analogue of the model. Suppose, that at the initial 
state we have an uniform distribution of energy of classical fields. The value of the energy density equals the absolute 
value of the static Casimir energy: po — n/2AL 2 . It corresponds to the classical potential A(t, x) — ip(t + x) + (p(t — x) 
with if(r) = 7tt/48L 2 . The classical energy is given by [631 ] : 



,L(t) n 

Edt) = J dxT m (t,x) = J 



t+L(t) 



t-L(t) 



(44) 



Next, we allow for the classical evolution of the electromagnetic system. From classical equations of motion we get: 
</?(t) = ¥?(/(t)). Using initial conditions, we obtain a classical global solution: 



tp(r) = 



48L 1 



R(t) . 



(45) 



We have encountered almost the same asymptotic formula for the energy as in the quantum case. The only exception 
is the coefficient: 



E c i(t) 



n tan 2 i 



12(2M- 1) 2 L 



(46) 



The results are confronted in Fig[6l We plot coefficients of asymptotic energy formulae for the first eight resonance 
channels. To make plots more readable, we have used continuous lines. Classically, the strongest growth is for the 
fundamental resonance frequency, next resonances are less effective. For the quantum case, the situation reverses. 
Due to the effect of quantum anomaly, there is no resonance in the first channel. Next, there is a rapid saturation 
for higher resonance channels. In units of 7rtan 2 9/12L, the sum of a classical coefficient and a quantum coefficient is 
just a unity. 
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4 5 

order of resonance 



Figure 6: Coefficients for quantum Eq. (|43|) and classical Eq. (|46[) asymptotic energy formulae. The coefficients are in units of 
7rtan 2 6i/12L. 



B. Homographic fundamental maps Ai 

We now turn the discussion to the case of maps Ai that are rational functions with single poles: 

h(v) 



A 1 (v) 



v-v 



(47) 



From Eq. (fTB|) we find that h(v) is at most a quadratic function. The periodicity condition Eq.© allows only for one 
singularity of Ai per period. Therefore, we can only consider homographic maps. It is convenient for us to start the 
discussion with inversions, and then we will be looking at a general case. 



Inversion map 



We evaluate that Eq. (fT7j) and other necessary conditions are satisfied by maps: 

Vq loL 



... Vn. UJ±j 

Ai(w) = , vq = tan-^- = tan I 



(48) 



There are no solutions for singular vo, so that u> ^ u>2N-i an d \9\ < ir/2. Moreover, we are forced to assume 9 ^ 
and this way all resonance frequencies are excluded here: to ^ lon- The auxiliary function / for r > L from Eq.0 
and its inverse function / _1 for r > —L yield: 

/(r) = ^arctan^/tan^f) + (|_JJ - 2M + 6(-0)) T , 
/ _1 (r) = £arctan(^/tan^) + ([^J + 2M + 0(0)) T , 

L=(M+f)T. 

The trajectory of the cavity wall for t > is reconstructed from Eq. (fT2"l) : 



(49) 



L(t) = L- 2 l + ^ 



arcsin (cos 29 cos {tot)) 



(50) 
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Evidently, the maximal velocity is now: 

Vmax = | COS26»| . (51) 

In the limit lj — > lun, we encounter a triangle wave trajectory. For u> = (lun +u>n-i)/2, our solution degenerates to a 
static one. The oscillations do always take place between L and L + sign#(7r — 4\0\)/uj. The corresponding milestones 
for our representation of Moore's function are given by: 

77 + 1 

L n = (-l) n L+ |_I-J(4M + signer . (52) 

We make the energy density explicit: 

~4lfe for r e [L 2 k-i,L 2k ) 

P(r) = { . „„..„. (53) 



uj 2 -uj v*(\+v 2 ) 

48tt 1 487T v^+v 2 



for t € [L 2k ,L 2k+1 ) 



There are wave packets in the energy density, but there is no unbounded growth of the total energy. The quantum 
cavity system is stable and its total accumulated energy oscillates with the period (4M + sign#)T. 

The solution is also well-defined for u < u>\ (L < T/2). It corresponds to M = and 9 > in Eq. p5|) and Eq. tfST))) . 
Here, as the only effect of the cavity motion there are pits in the energy density (negative wave packets) that may 
appear periodically in synchronization with cavity oscillations. 

2. Homographic map 

As a final and the most interesting application of our ideas, we consider a solution with a fundamental map being 
a homographic function. So then, upon confrontation with initial conditions Eq. (|17[) . we set: 

= _v 1 v + vo(vo~2v 1 ) 

V — Vi 

where vq — tan(wL/2) and v\ is an arbitrary parameter. It is straightforward to check that physical solutions exist 
on condition that Vq ^ v\. In passing, we note that the solutions that have been described in the previous subsection 
are reproduced for Vi = 0. 

The evaluation of relevant auxiliary functions / and / _1 ends with the results: 

/(r) = £arctaaAi(u) +(Lt + |J - 2M + Q{v-vi) - 0(«o - «i)) T , 
f-^r) = ^arctanAr^w) + + ±J + 2M + 0(« + Vj) - &{v! - v )) T , (55) 

L = (M + iarctani; )r . 

The milestones are given by: 



2 

L n = — arctan A n (v ) 



(2n + l)M + J2 e ( A fe l ( v o) + vi) - n0(«i - v ) 



fc=0 

The angle parameter 9 may be introduced here by using the following formula: 



T . (56) 



tan = i±^-7;o . (57) 
2v\ 

With the above definition, the derivation of the trajectory of the cavity wall from Eq. p^)) gives us: 

sin (ujL(t) + 9) = sin (uiL + 9) cos ut , (58) 

and it can be disentangled successfully: 

L(t) = L H [arcsin (sin (uiL + 9) cos (cut)) — arcsin (sin (uL + 9))] . (59) 
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We have assumed throughout this paper that the functions arcsin and arctan have their ranges restricted to 
[— 7r/2, tt/2). It makes the right hand side of Eq. ([55|) uniquely and properly defined. The maximal velocity is: 



\sm(ioL + 9)\ , 



(60) 



while the amplitude of oscillations is given by AL = (2/lu) arctan v max . 

It is important to point some special cases of Eq. ([59|) . For M = 1, Vq = and v\ = 1/(2 tan 8), we get a cavity model 
investigated by Law in [2l| . It was the first exact closed form solution presented in literature. The generalization of 
this solution for any M is leading to a second set of solutions described in the paper by Wu et al. [41| (the solutions 
with m = 2N in their notation). All solutions correspond to cavity vibrations with resonance frequencies u> — ujn- It 
was established for Law's and Wu's solutions that there appears resonant instability with a power-like behavior. We 
are not going to discuss these solutions here and send the reader back to the original papers. We wish only to note 
that there is one more class of solutions with a power-like instability. We obtain these solutions if we put vq — 2vi in 
Eq. (|54"]) . Let us describe them very briefly. The frequency of oscillation is a free parameter and may be tuned to any 
value greater than the fundamental frequency lo\ . But the amplitude of oscillations is uniquely defined by the choice 
of frequency. 

The solutions with exponential instability are obtained for v ^ and Vq =/= 2vi . The derivation of the maps A n for 
the representation of Moore's function Eq.Q requires the calculation of n-fold composition of the fundamental map 
given by Eq. ([5^|) . It is easy for homographies, so that the result is: 



A n (v) = A?» = (Ai - A 2 ) 



2(A' 1 ' + A 2 > + (A 1 ~A 2 )(A?~A 2 I ) 
4(AJ - \%)v + 2(A X - A 2 )(A» + A' 2 l ) 



(61) 



where: 



Ai, 2 = -ui ± vV)(2vi - v ) 
The profile function for the energy density is then given by: 



(62) 



P(r) = 



48tt 



48tt 



1 + v 2 



! ) 2 + («" + 



-) 2 



(63) 



We restrict ourselves to calculate only the approximate value of the total energy for large times. Therefore, the first 
term in Eq. (|63j) can be omitted, while the second term integrated over one period gives: 



dr p(r) 



uj\ Tr(H T H) 



48tt det(H) ' 

where H is a matrix composed of coefficients of the homography A„ in Eq. (f6"Tj) . It can be easily calculated that 



(64) 



Tr(H T H) _ 1 
det(H) 



4 V vi 



n(r) 



1 f 1 

2 V v i 



(65) 



For simplicity, we have assumed here that the number n(r) do not change in the interval of integration. Further, we 
assume small amplitudes: 2L(t) w 2L 2MT and obtain the approximate formula: 



E(t) 



AM 2 



1 



— ) cosli 

Vl 



y / V (2v 1 - Vq) t 

~ 2L 



Vl + ^Vq^V! - Vq) 



(66) 



In the above brief calculation, we have demonstrated that it is rather easy and safely in our treatment to perform 
approximate calculations and skip insignificant details. In fact, the treatment described in Section II and III is well 
adapted for perturbative methods. However, the relation for the minimal amplitude of oscillations AL m i n enough to 
trigger exponential instability of the cavity system vibrating at some fixed frequency uj can be derived exactly: 



AL r , 



UJ I _UJ_ _i_ 1 I 



(67) 



As the velocity of cavity wall cannot reach a speed of light, we obtain also the upper limit for amplitudes of oscillations 
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Figure 7: The phase diagram for stability and instability regions for the cavity model Eq. (|50|) . 



&-L min at given frequency u: 

AL 

max ^1 fnn\ 

r = — ■ (68) 

L to 

The above relation allows us to set up the phase diagram of stability and instability regions. A black area in Fig[7| 
correspond to the instability region of the vibrating cavity model Eq. (|50|) . Below this area, we have defined frequencies 
and amplitudes the cavity model is stable. A marginal behavior appears when the energy grows quadratically with 
time. It is observed for vq — (resonance frequencies, boundaries between adjacent Arnold's tongues) and Vq = 2v% 
(boundary points between stability and instability regions). Other parts of the diagram correspond to points where 
the cavity model is not well defined (physical assumptions about L(t) at the beginning of Section II are violated). 



V. CONCLUSIONS 



We have presented a rich class of exact and closed form analytical solutions for the quantum vacuum field in a 
one-dimensional cavity vibrating under the parametric resonance conditions. The solutions are valid for all times, 
frequencies of cavity oscillations and/or their amplitudes are free parameters. For small amplitudes, cavity oscillations 
are close to sinusoidal ones. In view of these properties, we can expect our solutions to yield all generic features known 
from other investigations on vibrating cavity models in a single dimension. 

The representation of solutions Eq.([5]) that appears in our treatment is based on SL(2, R) symmetry of scalar fields 
quantized in a static cavity. We have introduced the notion of fundamental maps that are more convenient to proceed 
than Moore's phase functions. There is a direct mathematical relationship between iterations of fundamental maps 
and the mechanism of parametric resonance. This is the way we can get insight into the regions of stability and 
instability of the model (see FigJTJ). One can calculate the rate of increase of the energy and the Lyapunov exponents. 
The stability-instability transition points and points between adjacent Arnold's tongues correspond to cavity models 
with a power-like instability. Thus, the most crucial questions can be tackled. If we insist on detailed calculations 
or exact formulas, then we have to set up how regions in space are covered by our maps (see Fig[T|). Summarizing 
technical aspects, we can tackle with any solution successfully and completely provided that we know its fundamental 
map Aj Eq.((7]) and respective ranges of maps L n together with n(r) Eq. ljlip . In this paper, general properties of 
fundamental maps for any physically reasonable solutions have been described. This setup is also a good start point 
for perturbative calculations. 
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To the best of our knowledge, both exact closed form solutions for off resonant frequencies and exact closed form 
solutions with exponential instability for vibrating cavities were not presented before. It refers also to Arnold's phase 
diagrams of stability-instability regions for solutions of vibrating cavity models. It is also important to state that the 
mechanism of parametric resonance in the quantum field theory shares many common features with its analogue in 
the classical theory. We believe that this similarity could be maintained also for three-dimensional vibrating cavities. 
The same should be true for the relevance of the symmetry of the quantized cavity model. 
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